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\D . Abstract 

The gravitational back-reaction on a certain type of rigidly-rotating cosmic 
^ ■ string loop, first discovered by Allen, Casper and Ottewill, is studied at the 

level of the weak-field approximation. The near-field metric perturbations 
"q | are calculated and used to construct the self- acceleration vector of the loop. 

Although the acceleration vector is divergent at the two kink points on the 
?h . loop, its net effect on the trajectory over a single oscillation period turns 

out to be finite. The net back-reaction on the loop over a single period is 
calculated using a method due to Quashnock and Spergel, and is shown to 
^ ■ induce a uniform shrinkage of the loop while preserving its original shape. 

The loop therefore evolves by self-similar evaporation. 

Short Title: Evaporation of a Cosmic String Loop 

PACS numbers: 04.25.Nx, 98.80.Cq 



1. Introduction 

Cosmic strings are thin filaments of Higgs field energy that may have played a 
role in the formation of large-scale structure in the early Universe. The dynam- 
ics of cosmic strings in Minkowski and Robertson- Walker backgrounds has been 
studied extensively [D El El IH El EJ , and their possible cosmogonic effects are well 



known (see [7j for a review), but their gravitational properties remain relatively 
unexplored. In particular, the motion of a string loop under the action of its own 
gravitational field is described by a non-linear system of equations that can nor- 
mally be solved only by numerical methods, even in the case of a zero-thickness 
string in the weak- field limit [B]. 

Quite apart from the potential cosmological ramifications, a better understand- 
ing of the gravitational field and back-reaction effects of general cosmic string 
configurations is important for the theoretical development of general relativity. 
Zero-thickness cosmic strings form a class of line singularities that are natural 
higher- dimensional analogues of black holes, and it remains an open question 
whether they can be satisfactorily incorporated into the framework of general rel- 
ativity. Areas that are usually regarded as problematic include the definition of 
an invariant characterisation of zero-thickness cosmic strings, the development of 
a suitable distributional formalism for line singularities, and the demonstration 
that zero-thickness cosmic strings are the unique limits of well-behaved finite- 
thickness solutions These issues can all be resolved in the few known cases 
of exact self-gravitating solutions [TUl HH H21 EE] , but a general resolution awaits 
further work. 

The uncertain status of line singularities in general relativity is underscored 
by the back-reaction problem. Almost nothing is known about the motion of 
a zero-thickness cosmic string under the action of its self-gravity. Although it 
can be shown that the self-gravity of a GUT string would almost everywhere be 
small enough to justify a weak-field treatment jH], even the weak-field dynamics is 
complicated and therefore little studied. Furthermore, flat-space string solutions 
commonly support pathological features - kinks and cusps - where the weak-field 
approximation is known to break down Much work needs to be done 

to justify the claim that zero-thickness cosmic strings can move freely and self- 
consistently under their own gravity, and are therefore compatible with general 
relativity. 

In this paper I calculate the gravitational field and back-reaction force on 
a certain type of rigidly-rotating cosmic string loop - the Allen-Casper-Ottewill 
solution |T5J - in the weak- field limit, and show that the loop evolves by self-similar 
shrinkage. This is the first explicit back-reaction calculation to be published 
for a cosmic string loop, as all previous calculations have relied on numerical 
approximation The importance of the results described below lies not so 
much in the simple behaviour of the loop, but in the detailed features of the 
gravitational effects involved. In particular, the loop supports a pair of kinks near 
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which the back-reaction force diverges but nonetheless remains integrable. 
2. Preliminaries 

A zero-thickness cosmic string is a line singularity which as it moves traces out a 
two-dimensional surface T whose dynamics and stress-energy content are governed 
by the Nambu-Goto action. In what follows, local coor- 

dinates on the four- dimensional background spacetime (M, g ab ) , and the metric 
tensor g ab has signature (+, — , — , — ), so that timelike vectors have positive norm. 
The world sheet T of a general zero-thickness cosmic string is then described para- 
metrically by a set of equations of the form x a = X a (( A ), where the parameters 
( A = (C^C 1 ) are commonly referred to as world-sheet or gauge coordinates. 

The intrinsic two-metric induced on T by a given choice of gauge coordinates 

is: 



material 7 AB is non-degenerate with signature (+, — ) almost everywhere on T, 
although ^ AB may be degenerate at isolated points called cusps. If 7 denotes 
I det^^g)! the Nambu-Goto action [TZl HHJ has the form 



where jj, is the mass per unit length of the string. 

The Lagrangian density £ = — /ry 1 / 2 in (12 .2)) is a functional of the position 
functions X a and their gauge derivatives, and of the background metric g a b. The 
equation of motion of the corresponding string is recovered by setting the func- 
tional derivative S£/SX a equal to zero, while the stress-energy tensor T ab at a 
general point on the world sheet T is ^g^^SC/Sgab, where g = \ det(g a fc)|. For 
the purposes of calculation, it is convenient to first fix the gauge so that 7 AB 
is everywhere proportional to the Minkowski tensor 77 AB =diag(l,— 1), a choice 
often referred to as the standard gauge. The corresponding gauge coordinates are 
conventionally written as C° = t arid C 1 = a. 

In the standard gauge, with X£ = dX a /da and X® = dX a /dr, the string 
equation of motion reads: 





(2.2) 



yti 1 -pa yo yc ya 1 pa 

A itt +1 fec A T A T - A , aa +1 bc A ff A^ 



(2.3) 
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where T^ c = ^g ad (gdb,c+9dc,b —Qbcd) is the Christoffel symbol on the background 
spacetime, evaluated at x a = X a (r, a). Also, the stress-energy tensor at a general 
spacetime point x c in the absence of other material sources is: 

T ab (x c ) = fxg- 1/2 J {X a T X b T -X a a Xl)5 {i \x c -X c )drda (2.4) 

The strong-field version of the back-reaction problem is completely specified by 
combining (|2.3|) and (|2.4|) with the Einstein equation G ab = —87iT ab , where G ab = 
R ab — \g ab R is the Einstein tensor. 1 

To date, the only known solutions to the strong-field back-reaction problem 
describe either infinite straight strings (with or without an envelope of cylindrical 
gravitational waves |T!H I2T] 122] ) or infinite strings interacting with plane- 
fronted gravitational waves (the so-called travelling- wave solutions (2S1I21])- One 
of the reasons why an explicit solution is possible in these cases is that the world 
sheets are intrinsically flat and non-dissipative. A self-gravitating string loop, by 
contrast, would be expected to continually radiate energy and angular momentum, 
and there seems little prospect of generating strong-field solutions of this type. 
The weak-field approximation is therefore the only tool currently available for 
analysing loop evolution. 

In the weak-field approximation g a b = r] ab +h a b, where r\ ab = diag (1, —1,-1, — 1) 
and the components of h a b are all small in absolute value compared to 1. To linear 
order in h ab the string equation of motion (|2.3|) then reads 

r, TT -r )m7 =Q n (2.5) 

where 

a a = (hl c -lv ad h bc , d )(X b X c a ~X b X c T ) (2.6) 

is the local four-acceleration of the string. Here, indices are everywhere raised and 
lowered using the Minkowski tensor i] ab and its inverse rf b . For future reference, 
h = h a a and inner products such as X%, X% and X a ■ X T are formed using r\ ab . 

If the harmonic gauge conditions h b a ,b = \h ia are imposed, the Einstein equa- 
tion G ab = —8nT ab reduces to nh ab = —16nS a b, where S a b = T ab — \r\ ab T^ and 

^^The convention adopted here for the Riemann tensor is that 

Rabcd = \ {gac,bd +gbd,ac -g a d,bc -gbcad ) + g e f {^aSbd ~ ^adFbc)- 

Also, geometrised units have been chosen, so that the gravitational constant G and the speed 
of light c are set to 1. 
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□ = d\ — V 2 is the flat-space d'Alembertian. With x a = [£,x] a the standard 
solution for h ab is: 

M*,x) = -4 / dV (2.7) 

y |x-x'i 

where t' — t — |x — x'| is the retarded time at the source point x'. 

When the source function S a b corresponding to the string stress-energy tensor 
(j2.4j) is inserted into (j2.7J) an integral over r, a and the three components of x' 
results. The distributional factor S^(x' c — X c ) can be integrated out by first 
transforming from y a = [r, x'] a to z a = [£', x'] a — X a (r, a), with a, t and x fixed. 
The Jacobian of this transformation is: 

\dy a /dz b \ = \dz h /dy a Y 1 = |X T ° - n' ■ X^" 1 , (2.8) 

where n' = (x — x')/|x — x'| is the unit vector from the source point x' to the field 
point x, and X a = [X°, X] a . On integrating over z a , t' is everywhere replaced 
by X° and x' by X. In particular, the gauge coordinate r becomes an implicit 
function of a through the light-cone condition 

X°(T,<r) = t-\x-X(T,(T)\. (2.9) 

The corresponding curve in the r — a plane, which represents the intersection of 
the past light-cone of the field point with the string world sheet, will be denoted 

by r. 

Since g 1 ^ 2 = 1 to leading order in h ab , the solution for h ab reads 

M*>x) = -4/i I 7^K°-n-X r |-W (2.10) 
Jr l x ~ x l 

where now n = (x — X)/|x — X| and 

^ ab (T, a) = X a T X b T -X a a X b a -\ V ab (X 2 - X 2 ). (2.11) 

The equation of motion (|2.5J) and the expression ()2.10j) for the metric perturba- 
tions h a b together constitute the weak-field back-reaction problem. 

It should be stressed that the weak-field equation of motion (|2.5|) holds only 
in a flat background, and in a particular gauge. Since questions have been raised 
about the validity of this equation [23] it is worthwhile taking the time to briefly re- 
view its derivation. The basic assumption underlying the weak-field back-reaction 
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formalism is that any solution pair (g ab ,X a ) to the strong-field problem repre- 
sented by equations (|2.3jl and (|2.4|) can - for small values of the mass per unit 
length n at least - be expanded as perturbative series of the form 

oo oo 

g ab = Y,H k g { $ and X a = Y,» k X( k ) (2-12) 

k=0 k=0 

where the functions and X? k s are independent of \i. 

In the problem at hand, g^ b = i] ab and figj^ = h ab , while the position vector 
X a appearing in ()2.5|) should be read as X^ + fiX?^ . However, for consistency the 
occurrences of X a in (|2.fij) . (|2.1()jl and (j2.11j) are all truncated at X? Q , only The 
weak-field equation of motion (|2.5J) is therefore effectively two equations, namely 

X(o),tt — -^(o))<7(7= (2-13) 

and 

5X a , TT -5X a , aa =a a (2.14) 

with 5X a = fiXfa. 

Equation ()2.5|) was generated by first fixing a gauge and then inserting the se- 
ries expansions (J2.12)) into the strong-field equation of motion (J2.3|) . An alternative 
method of derivation is to first write down the perturbative equation of motion in 
a gauge-independent form, and then fix the gauge. A suitable gauge-independent 
formalism has been developed by Battye and Carter , and in Appendix B it is 
shown that this second approach does indeed lead to the weak-field equations of 
motion (J2.13|) and (|2.14j) in the standard gauge. 

The structure of the paper from this point onwards is (I hope) very straight- 
forward. The unperturbed solution X^ is described in Section 0] (where, however, 
I will write X a in place of X^), the metric perturbation h ab is calculated in Sec- 
tion U?TJ the acceleration vector a a is calculated in Section |4~21 and the perturbed 
trajectory X a = X^ + 5X a is discussed in Section 14.31 However, many of the 
details of the calculation of a a are consigned to Appendix A. 

3. The Allen-Casper-Ottewill loop 

The standard method for solving the weak-field back-reaction problem, which 
was pioneered numerically by Quashnock and Spergel jS] over a decade ago, 
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is to start with a loop solution X a (r,a) to the flat-space equation of motion 
X a , TT —X a , acr = 0, calculate the corresponding self-acceleration vector a a at each 
point on the world sheet over a single period of oscillation of the loop, integrate the 
resulting weak-field equation of motion (|2.5J1 over a complete period to generate 
a new, perturbed solution of the flat-space equation of motion, and then iterate 
the procedure. What results is not a continuous solution to (|2.5jl but rather a 
sequence of flat-space solutions, each member of which is constructed by debiting 
from the preceding member of the sequence the gravitational energy and momen- 
tum it radiates over a single period. A quasi-stationary approximation scheme of 
this type is valuable whenever the mass per unit length \x of the string is small 
- such as is the case for GUT strings, which have // ~ 10~ 6 - as the sequence of 
flat-space solutions is then almost continuous. 

The Quashnock-Spergel method will be described in more detail in the next 
section. In this section I will briefly explain the choice of initial loop solution 
X a (r, a) to be used in the algorithm. The general form of loop solutions to the flat- 
space equation of motion X a , TT —X a , aa = is well known. In a flat background 
it is always possible to choose the gauge coordinates so that X°(r, a) = r. The 
spatial components of the position function X a have the form 



where for a loop in its centre-of-momentum frame the functions a and b are each 
periodic functions of their arguments with some parametric period L. The gauge 
constraints j TT + r ) aa = X 2 + X 2 = and j aT = X a ■ X T = together imply 
that | a' | 2 = |b'| 2 = 1, where a prime denotes differentiation with respect to the 
relevant argument. 

Since X(r + L/2, a + L/2) = X(r, a) for any values of r and a when a and 
b are periodic functions, the fundamental oscillation period t p of the loop is L/2 
rather than the parametric period L. The total four-momentum of the string loop 
on any surface of constant t = t is: 



p a = [i \ X^da = fi [1, ±a'(r + a) + |b'(r - a)] da = fiL[l, 0] (3.2) 



and in particular the energy of the loop (in this, the centre-of-momentum frame) 
is E = [iL. The corresponding angular momentum of the loop is: 








(3.3) 
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In the weak-field approximation the gravitational power per unit solid angle 
radiated by a periodic source out to a point x in the wave zone, where r = |x| is 
much greater than the characteristic dimension L of the source, is given by 

% = ~£t (3-4) 

where n = x/r and R ab here denotes the Ricci tensor evaluated to second order 
in the metric perturbations h a b and to second order in r _1 . The angled brack- 
ets indicate coarse-grained averaging, which is performed by expanding h a b as a 
Fourier series in the components of the field point x a and then eliminating all but 
the zeroth-order Fourier mode from R ab . 

The resulting expression for the power radiated per unit solid angle is 

rIP i , 2 00 

^ = 7E m2 l r tr^ 2 l (3-5) 

m=l 

where 10 = 2it/t p is the circular frequency of the source, and 

f a \m,Yv) =t~ l [ I" T ab {t,x)e imu}it - nx) dtd 3 x. (3.6) 

</R 3 JO 

In the case of the string stress-energy tensor ((231), with xa = l T > |( a + b )] a 
a given flat-space solution, it is more convenient to work with the light-cone 
coordinates a± = r ± a rather than with r and a themselves. The fundamental 
domain < r < L/2 and < a < L is then equivalent to < a + < L and 
< cr_ < L. If a' c = [l,a'(a + )] c and b' c = [l,b'(a_)] c then f ab = fiLA^ a B b \ 
where the round brackets denote symmetrisation, and 

A c (m,n) = ±[ e 27Tim[a +- n - a( - a+)]/L a' c (a + )da + (3.7) 
Jo 

and 

5 c (m,n) = i/ e 2nim[a — n - h{a - )]/L b' c (a_)da_. (3.8) 
Jo 

In terms of the four-vectors A c and B c , therefore, 
dP 00 

= Stt/x 2 m2 i( A ■ A *)( B ■ B *) + \A-B*\ 2 -\A- B\ 2 }. (3.9) 

m=l 
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A convenient measure of the total power radiated by a string loop is the 
radiative efficiency 7 = [i~ 2 P, where P = f ^dVt is the total power of the 
loop. In terms of 7 , the fractional energy radiated by the loop over the double 
oscillation period L = 2t p is \ AE/E\ = (PL)/(pL) = /Z7 . Numerical simulations 
of networks of string loops indicate that 7 is typically of order 65-70, although 
it may fall as low as 40 [27, 28] • 

The range of 7 has been examined analytically by Allen, Casper and Ottewill 
16 for a large class of string loops in which the mode function a consists of two 
anti-parallel segments of length aligned along the z-axis, so that 

a (^) = ( (L-a + )z if i L <a + <L> (3 ' 10) 

and the mode function b is restricted to the x-y plane but is otherwise uncon- 
strained. The minimum possible value of 7 for loops of this class is approximately 
39.0025, and occurs when 

b(tr_) = ^[cos(2ttct_/L) x + sin(27ra_/L) y]. (3.11) 

The corresponding loop consists of two helical segments which are rigidly ro- 
tating about the z-axis with an angular speed to = 4n/L. The evolution of this 
loop is illustrated in Figure 1, which shows the y-z projection of the loop at times 
r-e = 0, L/16, L/8 and 3L/16 (top row) and r-e = L/4, 5L/16, 3L/8 and 7L/16 
(bottom row), where the time offset e is 0.02L. The string has been artificially 
thickened for the sake of visibility, and the z-axis is also shown. The projections 
of the loop onto the x-y plane are circles of radius L/(4ir). The points at the 
extreme top and bottom of the loop, where the helical segments meet and the 
modal tangent vector a' is discontinuous, are technically known as kinks. They 
trace out circles in the planes z = L/4 and z — 0. All other points on the loop 
trace out identical circles, although with varying phase lags. The net angular 
momentum of the loop is J = ^fiL 2 z. 

The importance of the solution defined by the mode functions ()3.10|) and 
(|3.11|) . which I will henceforth call the Allen-Casper-Ottewill (ACO) loop, lies in 
the fact that its radiative efficiency 7 ~ 39.0025 is the lowest known of any loop 
solution. Casper and Allen |2Hj have examined an ensemble of 11,625 string loops, 
generated by numerically evolving a large set of parent loops forward in time (in 
the absence of radiative effects) until only non-intersecting child loops remained, 
and found that only 6 of the loops studied had radiative efficiencies less than 42, 



9 



while none at all had 7 < 40. When examined more closely, all 6 of the low-7 
loops were seen to have the same general shape as the ACO loop. A second study 
by Casper and Allen involving another 12,830 loops yielded similar results. 

This observation has interesting implications for the evolution of a network 
of string loops. If it can be shown that an ACO loop does not evolve secularly 
towards another shape with higher radiative efficiency, then it will be the longest- 
lived loop in any ensemble of loops with a given energy E (and invariant length L). 
Furthermore, if the general trend in loop evolution is from higher to lower radiative 
efficiencies - a trend which has not so far been demonstrated, but which would be 
expected if radiative dampening acts to preferentially eliminate higher-multipole 
structure - then the ACO loop would in some sense play the role of an attractor in 
the ensemble of loops forming a network. However, it is conceivable that, because 
the energy E of a string loop falls off as Eq — fi 2/ ~f°t (since |Ai?| = fi^E over a 
time period At = 2t p , but At = E/fi), any fall in the radiative efficiencies of the 
loops in a network might proceed too slowly to have a significant effect on any 
but the very largest loops. 

What I will demonstrate below is that, in the weak-field approximation at least, 
the ACO loop does indeed evaporate in a self-similar manner, and so 7 remains 
constant throughout its evolution. When the weak-field equation of motion (|2.5|) 
has been integrated over a single period t p , it turns out that an initial ACO loop 
with energy E and invariant length L decays to form a loop with exactly the same 
shape, but with smaller values E' = E(l — h(J/y°) and L' = L(l — 5^7°) of energy 
and length, plus a small rotational phase shift. 

4. Back-reaction of the ACO loop 
4.1. The metric perturbations 

The first step in the back-reaction calculation is to evaluate the line integrals (|2.7|) 
for the metric perturbations h^. Let X a = [r, X] a be a general source point on 
the string loop, and x a = [i, x] a a general field point. The mode functions of an 
ACO loop with invariant length L are as given in ()3.1U|) and ([3.11jh but can be 
simplified slightly (by rescaling the gauge coordinates a + and a J) to read 



a(cr + ) = L\v + \ z if — \ < v + < \ 



(4.1) 



and 



b(cr_) = ^[cos(27rt>_) x + sin(27rt>_) y]. 



(4.2) 
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where (v + ,V-) = (er + /L, er_/L). Note that at any point on the world sheet it is 
always possible to choose v + and t>_ so that r = \L{v + + u_) remains equal to 
the coordinate time t, but v + G [— 4, |]. 

The spacetime coordinates (|3.1j) of the source point are then 

X a = \L[v + + u_, ^ cos(2tto_) x + ^ sin(27ro_) y + |v + | z] a . (4.3) 

For a general field point x a , with x = xit + yy + z 2, the light-cone condition 
i — r = |x — X| reads: 

2t/L-v + -V- = {[2x/L-± cos(2vr ? ;_)] 2 +[2y/L-^ sm(2nv_)} 2 +(2z/ L-\v + \) 2 } l l 2 , 

(4.4) 

or equivalently 

v+ = -§(v_ - 2i/L - 2sz/L) 

m\2x/L - i cos(2tto_)] 2 + [2j//L - ^ sin(2vn;_)] 2 }/(t;_ - 2t/L + 2sz/L) 

(4.5) 

where s = sgn(t> + ). 

In particular, if the field point itself lies on the loop, and has scale-free gauge 
coordinates (a + /L, &-/L) = (u+,U-) so that 

[i, x] = \L[u+ + U-, cos(27ra_) x + ^- sin(27ra_) y+|u + | z] (4.6) 

then the light-cone condition (|4.4j) reads: 

u + - v + + m_ - v_ = [^{1 - cos2tt(u_ - v_)} + (fu+ - sv + ) 2 } 1/2 (4.7) 

where s = sgn(w + ). [For future convenience, the 4- vector on the right of ()4.6j) will 
abbreviated as X a = [r, X]".] 

Recall from the previous section that the ACO loop consists of two helical 
segments bounded by kinks at a+/L = — 5, and |. If the source point and field 
point lie on the same segment (ss = 1) then su + — sv + = ±(w + — v+) and (|4.7j) 
can be rearranged to read 

2(u_ - - u+) + (u_ - vJ) 2 = -^{1 - cos2tt(m_ - «_)}. (4.8) 

For a fixed field point (w + ,-u_) this equation clearly has two possible solutions: 

w_ = u_ (4.9) 
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and 



v + - u + = - u_) + ^{1 - cos27r(t;_ - u_)}/(i;_ - u_). (4.10) 

In view of (|4.7|) the first solution is feasible if and only if v + <u + . It is also easily 
verified that if the second solution holds then u + — v + + U- — v _ > if and only 
if v- — u- < 0, in which case ()4.10|) indicates that v + > u + . So the two solutions 
are mutually exclusive and fully cover the case where the source point and field 
point lie on the same segment. 

If on the other hand the two points lie on different segments, then ss = — 1 
and su + — sv + = ±(w+ + v + ). So (|4.7j) expands to give: 

(2u + - t>_ +U-)(-2v + - v- = 2^r{l - cos2tt(w_ - v-)}, (4.11) 

or equivalently 

v+ = - u_) + ^{1 - cos27r(w_ - m_)}/(w_ - u_ - 2u+). (4.12) 

For fixed values of w + and m_ this equation typically defines two branches in the 
v + -v - plane, one with V- — w_ — 2m + < and the other with t>_ — w_ — 2m + > 0. 
From (|4.1ip . the combinations 2u + — v _ + m_ and — 2f + — u _ + m_ are either 
both positive or both negative, and since their sum 2{u + — v + +u_ —v_) must be 
positive or zero from (|4.7jl . the branch with f _ — u_ — 2m + < is the only feasible 
one. 

Figure 2 shows the backwards light cones of two field points on the same 
spacelike section f = of the loop. The field point marked A coincides with the 
kink at u + = u_ = 0, and its past light cone consists of the straight line segment 
f_ = (solution (|4.9Jl ) on the left, and the curve v + = — ^ v - + ^2 (1 — cos27rw_)/w_ 
(which can interchangeably be regarded as solution (|4.1(J|) or solution (J4.12j0 on 
the right. The point marked B is a general field point, and its past light cone 
contains all three components identified above: solution ()4.12|) on the left, followed 
by the straight-line segment (|4.9)1 . and then solution 1)4.10)1 on the right. For future 
reference these three segments will be designated T~, T and T + respectively. 

It is evident from the foregoing discussion that the integration contour T is a 
piecewise-smooth spacelike or null curve that wraps once around the world sheet 
T of the loop. When calculating the metric perturbation or its spacetime 
derivatives h ab , c at a given field point [i, x] it is possible in principle to use a or v + 
or f_ as the integration variable on T, as circumstances warrant. From ()2.9j) and 
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the fact that a = \L(v + — vJ) and t — r = |x — X|, the relevant transformation 
rules are: 

|x - X| -1 |l - n ■ X T | -1 d<r = LAZ 1 dv+ = -LA^dv^ (4.13) 

where 

A± = t-r-2(x-X)-X± (4.14) 

and it is understood that v + is varied with t>_ held constant, and vice versa. 
Because an explicit formula (J4.5|) exists for v + as a function of t>_ on T, the most 
useful choice of integration variable will normally be t>_, in which case (I2.1U|) 
becomes 

h ab (i,5t) =AfxL J A-^^dv^ (4.15) 

Now, in terms of the derivatives of the light-cone coordinates a + and cr_ the 
formula (|2.11|) for the source function ^> ab reads: 



ah 



2(X1 X"_ + X a _ X b + ) -2 7] ab X + -X 



(4.16) 



where X± = [|, X ± ] a , and X + = |s z and X_ = |(— sin 2nv _ x + cos 27ru_ y) are 
orthogonal vectors. Hence, X + ■ X_ = | and 



Also, since 



1 

sin 2ttv_ 

— COS 27TV- 
— S 



sin 27ru_ 
1 


-s sin 2ttv _ 



— COS 27TU_ 


1 

S COS 27Tf _ 



- LA; 1 = 2(w_ - 2t/L + 2SZ/L)- 1 
the integral ()4.15|) for h ab reads 



—s sin 27rv_ 

S COS 27TU_ 

1 



(i 6 



2t/L + 2sz/L)- i ^ afe cfe 



(4.17) 
(4.18) 

(4.19) 



where the limits on u_ are fixed by the requirement that v + range from — ~ to 
(with s = — 1) and then from to ~ (with s = 1). Note here that in view of (J4.13|) 
and the fact that A + and A_ are non-negative functions, v _ is a non-increasing 
function of v. on T. 
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If V-i, V and V x denote the values of t>„ when v + = — |, and | respectively, 
then the equation (|4.4|) for T indicates that 

2t/L-V k -l = {[2x/L-^cos(27rV k )} 2 + [2y/L-^ S m(2nV k )] 2 + (2z/L-\l\) 2 } 1/2 

(4.20) 

in each of the three cases k = —1, and 1. In particular, it is evident that if V x 
satisfies P~2"Uj) with k = 1 then V_i = V 1 + 1 satisfies flQUP when fc = -1. So 
Vi < V < V-i = V x + 1, and (ETToT) becomes 

/•Vo rVi+1 
h ab (t, x) = 8/i / (t>_ - ^_)~ 1 ^'afe|s=i dv- + 8/i / (f _ - ^ + ) _1 ^ afe | s= _i dv^ 

JVi JVo 

(4.21) 

where ip± = 2t/L ± 2z/L. 

In the case of greatest interest, where the field point [t, 5c] is a point (j4.6|) on 
the loop, if u + > then Vo = w_ and (from (j4.7|0 

V\ = M _ + W x (u + ) where u+ = ± + \W X - ^(1 - cos 27rWi)/Wi (4.22) 

while if u + < then V x — it_ — 1 and 

Vo = u_ + Wo(u+) where w+ = §Wo - 4^2 (1 - cos 2ttWo) /W . (4.23) 

Note that W x G (-1, 0) for u+ G (0, §) and W G (-1, 0) for u + G (-|, 0). In fact, 
W x (u + ) = Wq(u + — ~) for all u + G (0, |). At each of the kink points u + = — ~, 
and I the values of (VI, Vo) tend to (w_, m_) as u + approaches the kink value from 
below, and tend to (w_ — 1, uJ) as u + approaches the kink value from above. 
Collecting together the results of the preceding paragraphs gives the formula 

rV6— 4>- rVi+l-4> + 

h ab (i,x) = 8/i / U7 _1 ^ ab (l, w+4>_) dw+8/j, / w' 1 ^ ab (-l,w+ip + ) dw 

Jv 1 -i>_ Jv -ip + 

(4.24) 

where the components (j4.17J) of ^ ab are all smooth, and all except ^ xy and ty yx 
are non-zero. The two integrals in ()4.24|) are guaranteed to converge if the field 
point [i, x] lies off the string, but not if it lies on the string. In the second case, 
ip± = u_ + u + ± \u + \ and so V — ip_ = if u + > while V x + 1 — ip + — if 
u + < 0. Thus at least one of the integrals in ()4.24|) diverges. 

To investigate more fully the behaviour of the metric perturbations near the 
world sheet T, suppose that 

[t,x] = [r,X] + [Si, 65c) (4.25) 
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where 8t and <5x are assumed to be small. Then if u + > the equation (|4.2U|) 
with k = can be solved for Vq to second order in [5i, 5x] to give 

Vo-tf)_ w -^2{5s 2 + 5y 2 + (5t-5z) 2 + 2(5t-5z)[fesin(27r?i_)-5ycos(27TM_)]} 

(4.26) 

whereas if it + < the equation (j4.20j) with k — 1 can be solved for Vi to give 

\4+l--0, « - 1 2 L, T2 {5x 2 +5y 2 +(5i+5z) 2 +2(5i+5z)[5x sm(27ru_)-5y cos(2tto_)]}. 

(4.27) 

Here, the quadratic forms 

Q ± = Sx 2 + 5y 2 + (Si ± 5z) 2 + 2(<Jf ± 5z) [5x sin(27rn_) - 5y cos(2ttu_)] 

= 2^ ab (Tl, u^[5i,5R} a [5i,5R] b (4.28) 

are positive semi-definite and vanish on the 2-surface in [5i, 55c] space spanned by 
the vectors [1, =Fz] and [1, — sin(27ru_)x + cos(27rw_)y], which is just the tangent 
plane to the string at the point [f, X]. Thus, if the field point [i, x] approaches 
a point on the string from any direction outside the tangent plane, the divergent 
contribution to h a b has the form 

/i a6 ^8/i* a6 (±l,^)lnQ T (4.29) 

where the upper sign applies when u + > and the lower sign when u + < 0. Note 
also from (J4.26j) and (J4.27|) that the divergence is more complicated at the kink 
points u + = and u + = — |. 

The divergence of the metric perturbation h a b on the world sheet is to be 
expected (see for example jjj, p. 216), and is simply a reflection of the conical 
singularity that is present at every ordinary point on a zero-thickness cosmic 
string. For example, the exact metric about a static straight string lying along 
the z-axis has the "isotropic" form j2Hl 

ds 2 = dt 2 - dz 2 - p^(dx 2 + dy 2 ) (4.30) 

where p = (x 2 + y 2 ) 1 ^ 2 . In the limit of small fx the only non-zero metric per- 
turbations are h xx = h yy = 8/ilnp, which satisfy the Lorentz gauge conditions 
h b a ,b = \h, a but also diverge logarithmically as p — > 0. 

Two quantities that will play an important role in Section 14.31 in fixing the 
initial data for the perturbation 5X a are the values of the functions habX+X 1 ^ 
and KbX^Xt on the world sheet T. Since ^ a bX b + = ^ a bX b _ = at all ordinary 
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points on T, the leading-order divergence (J4.29|) in h a b makes no contribution to 
these functions. 

At the next highest order in [5i, 5x] , the metric perturbation has the form 



hab = 4/i 
where, for < u + < \ 







s + + 




-C+ - C- 


-I+- 


f/_ 




-S- 


I+ + 


J_ 





-s+- 


fS- 


-c+ 


-c_ 











C_ 


-h 


+ 1- 


-s+- 


-s. 


. c+ - c_ 


I+ + 


J_ 


< I 

- ~ 2' 














h 




lnl^l, 




J_ = In D+ - 


In 2u + 





J aft 



(4.31] 



(4.32) 



(5 +J C+) = (cos27rw_,-sm2™_)Si(27r|Wi|) 

- (sin 2ttm_ , cos 2ttm_) [Ci(2vr |Wi|) - 7e - ln27r] (4.33) 

and 

(5L,C_) = (-cos2tt(2m + + M_),sin27r(2M + + M_))[Si(27rD + ) - Si(47TM+)] 
+ (sin2vr(2M + + m_),cos2vt(2m + + M_))[Ci(27rL) + ) - Ci(4vra + )] 

(4.34) 

with Si(x) = J* w -1 sin wdw and Ci(x) = — J x °° w^ 1 cos w dw the usual sine and 
cosine integrals, and 7 E ~ 0.5772 Euler's constant. (In fact, Ci(27ra;)— 7 E — In 2ir = 
lnx — Jq x~ x (1 — cos 2irx)dx.) Also, the function D + {u + ) is defined by 

2itD + = 2n(2u + - Wt - 1) = -J(l - cos27rWi)/Wi (4.35) 

while Wj is defined implicitly as a function of w + through (|4.22j) . 

Substituting X" = |[l,sz] a and X" = |[1, — sin 2%U- x + cos 2-ku- y] a gives 

h a6 X£X + = 4/i/_ = 4/i(ln D+ - In 2u + ) (4.36) 

and 

h ab X a _X b _ = 2/i(/+ + J_) - 2fi[(S + + SL) sin 2ttu_ + (C+ + C_) cos 2vra_] 
= 2/i(lnL> + - ln2u+ - ln|Wi|) + 2//[Ci(27r|Wi|) - 7 B - m27r ] 
-2//[Si(27rD + ) - Si(4vra + )] sin(4vra + ) 

-2fi[Ci(2nD + ) - Ci(47rn + )] cos(4t™ + ) (4.37) 
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when < u + < \. The function h a bX^_X^ diverges as 4/ilnw + as u + — > 0, and 
as g/iln(l — 2u + ) as u + — > \. However, the function h a bX"iX b _ is bounded, and 
converges to 2/x[Ci(27r) — 7 E — ln27r] « —4. 875 3/x at both w + = and m + = |. 

In the case where — | < u + < 0, it turns out that J + (« + ) = /_(«+ + |) and 
I-{u + ) = I+(u + + |), and analogous identities hold for S± and C±. The functions 
h a bX+X+ and h a bXtX b _ are therefore periodic functions of w + with period |. 

A similar calculation of the leading-order contributions to h a b near the kink 
point u + = gives: 

5+, C_|_) ~ (1, sin 27rti_, cos 27tm_) ln(L _1 Q-/(5o) 
+ (0, cos 27tm_, — sin27TM_) Si(27r) 

-(0,sin27ru_,cos27TM_)[Ci(27r) -7 E -ln2vr] (4.38) 

and 

(I-,S-,CJ) w (1, sin 2ttm_, cos 2ttu_) ln(4L _1 Q ) (4-39) 

where 

Q = Si + 5x sin 2itU- - 5y cos 2tt?i_ = 2r] a6 X!; 5x] 6 . (4.40) 

Hence, habX^X_ = 2yu[Ci(27r) — 7 E — ln27r] at the kink point and the function 
habXtX'L is continuous there. By contrast, X" has a jump discontinuity at the 
kink point, and h a bX+X+ (as well as containing formal logarithmic divergences) 
is undefined even for non-zero [Si, 5x] . 

4.2. The self-acceleration vector 

Although h a b diverges at all points on the loop, the self-acceleration vector a a 
defined in (12. 6|) typically does not. This result was first established by Quashnock 
and Spergel for a generic point on a smooth loop, but does not extend to non- 
generic points such as cusps and kinks. 2 In terms of the light-cone derivatives X" 

2 It should be emphasised that the Quashnock-Spergel result holds only to linear order in 
the expansion parameter /i. Copeland, Haws and Hindmarsh |30| have made the claim that 
the integral for the self-acceleration vector a a contains a local divergence of the form /z 2 In p, 
where p is the normal distance from the world sheet. Carter and Battye |25j have more recently 
disputed this claim, maintaining that the divergence is cancelled by a counter-term that was 
previously neglected. Since neither analysis is based on a full second-order expansion of the 
equations of motion f l|2.3() and the equation K c = of Appendix B in the respective cases), this 
dispute remains unresolved. However, the question has no relevance for the purely first-order 
treatment of the back-reaction problem offered here. 
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and Xf. of X a = [r, X] a , the expression (J2.6j) for a a at a point on the 

string loop can be rewritten as 



a>+,ti_) = -r j ad (2h db , c -h bc , d )(X b + X c _ + X*_X_ 



(4.41) 



where it is understood that the right-hand side is evaluated by taking the limit 
[i , x] — > [f , X] , from any direction outside the tangent plane, in the metric deriva- 
tives h ab , c . 

Now the ACO loop, as described by ()4.1|) and ()4.2j) . rotates uniformly about 
the z-axis, and so the self-acceleration vector a a is uniformly rotating as well. In 
order to simplify the calculation of a a , it is useful to premultiply a a by the matrix 



10 

cos27ru_ sin27ru_ 

— sin 27ru_ cos 2nu- 

1 



(4.42) 



ab 



the effect of which is to map each horizontal radius vector joining the rotation 
axis to a point on the loop into a vector parallel to the x-axis. In other words, 
a L = M£a b is the lateral component of the loop's horizontal acceleration and 
a N = M^a b is its normal component. 

To demonstrate that a a is well defined away from the kink points at u + = 
and ±|, and to explicitly calculate the components of the co- rotating acceleration 
vector M^a b , is straightforward but tedious. The details can be found in Appendix 
A. It turns out that a 1 , a L and a N are periodic functions of u + , and a z is an anti- 
periodic function, with period ~. 

For < u + < |, the components of M^a b are 



a 



a 



32TT 2 fiL- 1 F(W 1 )W 1 (l - cos 2irW x ) 



(4.43) 



and 



a 
a 



N 



—32irfiL' 



{Si(47ra + ) - Si(27rD + )}sin47rn + + {Ci(47ra + ) - Ci(27rD + )} cos4t™ + 
-{Si(47TM + ) - Si(27rD + )}cos4vTM + + {Ci(47rw + ) - Ci(27rD + )} sin47™ 

2ttW 1 cos 2irW 1 - sin 2irW 1 
cos 2-kW! - 1 + 27rW! sin 2^! 



+327r 2 fxL- 1 F(W 1 )W 1 



(4.44) 
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where 

F(Wi) = (2n 2 W? + 1 - cos 2ttWi - 2nW x sin 2nW x y x , (4.45) 

27TD+ = -i(l - cos27rWi)/Wi. (4.46) 
as before, and W\ is defined implicitly as a function of u + through 

Airu+ = 2tt(1 + Wi) - J(l - cos27riy 1 )/W / 1 (4.47) 

(see 

The dependence of the acceleration component a 1 on u + is shown in Figure 
3. The interval u + G (0, |) corresponds to Wi G (—1,0). Since a* diverges as 
32/iL" 1 W / 1 " 1 near Wj. = while u + — ~ « I^Wf for small values of VFi < 0, the 
singularity in a* at w + = has a 1 ps — 32(|7r 2 ) 1 / 3 / uL _1 /|'u + | 1 / 3 for w + < 0. Note 
however that the singularity in a 1 is one-sided, as a 1 ~ — 327r 2 / uL~ 1 (l + W^i) 2 for 
W 7 ! near — 1, while u + « ~(1 + Wi), and so a* w -^STr 2 /!^ -1 ^^ for n + > 0. 
Similar remarks apply to a z , save that a z > for u + G (— 1,0) and a 2 < for 

w+e(o,|). 

The point m + = r + a = of course marks the location of the lower of the 
two kinks on the string, which propagates around the string in the direction of 
increasing a. For any fixed value of f , points with u + > lie ahead of the kink and 
have yet to be disturbed by its passing, whereas points with u + < lie just behind 
the kink and so have recently been subject to severe gravitational stresses. This 
is the reason for the one-sided nature of the singularity in a 1 and a z at u + = 0. 
Also, a 1 < because the string is radiating energy, whereas a z has different signs 
on the two branches of the loop because it is directed counter to the ^-component 
of the string's 4-velocity X", which is negative for u + G (— 1,0) and positive for 

u + e {().,{). 

The two remaining components of the acceleration vector, a L and a N , are 
plotted in Figures 4 and 5 respectively. Both are divergent at u + = 0, with 
a L ~ y^/xL -1 In for u + < and a L ~ 32iTfiL~ 1 lnu + for u + > 0, while 
a N w -32(|tt 2 ) 1 / 3 / uL-7|m + | 1 / 3 for u + < and a N w 1287r 2 /iL- 1 M+ lnw+ for 
u + > 0. The normal component is negative almost everywhere because it 
acts to torque down the rotational motion of the loop, and like a 1 and a z it is 
strongest in the immediate wake of the two kinks. The lateral component a L 
is also negative, but its singularity at u + = is more symmetric. The lateral 
acceleration is presumably a response to the imbalance between the string tension 
and the centripetal force caused by the rotational torquing. Its net effect is to 
induce lateral shrinkage of the loop as it radiates. 
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Although all four components of the acceleration vector a a are divergent at 
the kink points u + = — |, and ^, the singularities are all integrable, in the sense 
that J a a du + exists on any closed sub- interval of [— |, |]. This is turn means that 
the weak-field equations of motion for the ACO loop are integrable as well, as will 
be seen shortly. 

It could be argued that the method adopted here, of evaluating a a away from 
the kink points then integrating the equations of motion through the kinks, ig- 
nores the possible presence of other localised (delta-function or non-integrable) 
singularities in a a at the kink points themselves. There seems to be no way of 
ruling out this possibility short of solving the full set of field equations for a finite- 
thickness string and taking the limit as the thickness goes to zero. However, it 
should be noted that a formal calculation of a a at the kink point u + = us- 
ing only the logarithmically divergent contributions (|4.38j) and (|4.39jl to h ab gives 
a a = irrespective of the value chosen for s in the tangent vector = |[1, sz] a . 
(This remark should not be taken to mean that a a actually vanishes at the kink 
points, as terms of order [5i, 85c] in h ab contribute additional terms to a a which 
have no well-defined limits when [8t, 55c] —>■ 0. 3 ) Furthermore, as will be seen, the 
energy loss of the loop over a parametric period exactly matches the radiative flux 
calculated for the ACO loop in JH], so there are good reasons for believing that 
the acceleration vector a a has been represented accurately. 

4.3. Secular evolution of the ACO loop 

Once the acceleration vector 

a a = [a*, (a L cos 27ra„ — a N sin 27ra_)x + (a L sin 2nu_ + a N cos 27ra_)y + a z z] a 

(4.48) 

3 More explicitly, an expansion of h a b about the kink point u + = reads, to linear order in 
8x a = [Si,Sx} a , 

h ab = 8M[* a6 (l,u_)ln(L- 1 Q_/Oo) + *a b (-l,u-)ln(4L- 1 Q ) + ^afc(w-) 
+F- b (u-){2(St/L - Sz/L) ln(L- l Q_/Q ) - L^Q^/Qo} 
+F+(u^){2(Si/L + Sz/L) HAL^Qo) + L~ 1 Q + /Q Q ] + • • •] 

where the ellipsis (• • •) includes terms of first-order in 8x a even more complicated than those 
shown. The corresponding acceleration vector a a , calculated from l)4.41[l . contains no contri- 
butions from the first two, logarithmically divergent, terms. However, the terms proportional 
to F~ b and Fi, and the other linear terms, do formally contribute to a a , with coefficients that 
typically depend on s = sgn(u + ). The derivatives with respect to Sx a of all these terms are 
undefined in the limit as Sx a — > 0. 
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is known as a function of u + and it_, it is a straightforward matter to use (|2.5|) 
to calculate the secular perturbations in the trajectory of the ACO loop due to 
gravitational back-reaction. In terms of the derivatives with respect to the scaled 
gauge coordinates cr + = Lu + and cr„ = Lu_ the first-order equation of motion 
(gZEB reads 

6X, a + _=la a (a + ,a^). (4.49) 

The initial data for this equation on the curve 4 r = |(cr + +o~_) = are fixed by 
requiring £X+| T=0 and 5X^\ T=0 to satisfy the first-order gauge constraint (|B.13|) . 
while preserving the initial configuration X a | T= o of the loop and the tangent space 
to the loop at each point on the initial curve. This can be done by writing 

6Xl\ T=0 = k + Xl\ T=0 + l + X% =0 (4.50) 

and 

5X a _\ T=0 = k~X a _\ T=Q + rXl\ T=0 (4.51) 

where the coefficients k^ and ^ are functions of a = |(cr + — cr_) to be determined, 
while 

X^| r=0 ( ( r) = |[l,sgn(a)z] a (4.52) 

and 

X% =0 (a) = i[l,sin(27T(j/L)x + cos(27r(T/L)y] a . (4.53) 
In the light-cone gauge, the constraint (jB.13|) reads simply 

h ab X a + X\ + 27] ab 5XlX b + = and h ab X a _X b _ + 2r] ab 5X^X b _ = (4.54) 

and (since rj^XtX^ = \) is satisfied on the initial curve r = if 

l + = -2(h ab XlX b + )\ T=0 and V = -2(h ab X a _X b _)\ T=0 (4.55) 

where hahX+X 1 ^ and habX^X^ are known functions of u + = (r + a)/ L given by 
fl4.3fij) and ()4.37|1 respectively. 

The initial configuration of the loop will be preserved if 

= 5XX =0 = 5X a + \ T= v - 5X% = v = (k + - r)X»| T=0 + (/+ - AOJT | T=0 (4.56) 

at all points on the initial curve, and so k + = l~ and / + = k~ . Together with 
8X a \ T= Q = 0, the initial data for the perturbation 5X a therefore have the form 

5X a + \ T=0 (a) = 5X a _\ T=0 (a) = -2(h cd X c :X d )\ T=0 Xl\ T=0 - 2(h cd X c + X d + )\ T=0 X% =0 . 

(4-57) 

4 The bars on X a , f and a can henceforth be dropped without risk of confusion. 
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It might seem at first glance that replacing X±\ T=0 with X±\ T=0 + 5X±\ T=0 
will induce a change in the initial velocity of the loop. However, the effect of the 
changed initial conditions is to simply realign the gauge coordinates in response 
to the rotation of the light cones caused by the metric perturbation h a b- In the 
standard gauge, the tangent vector X% at any point P on a given spacelike section 
r = constant of the world sheet is uniquely determined by the tangent space at 
P and the choice of gauge coordinate a on the section, as X" is just the future- 
pointing projection of the tangent space orthogonal to X%, normalised so that 
X 2 = -X 2 a . 

The equation of motion f|4.49j) can now be integrated forward to give 

5X a + (a + ,a_) = \ ^ a a (cr +J 6)de + 5X a + \ T=0 (a + ) (4.58) 



and 



5X a _(a + ,a_) = \ I + a a (6,a„) d6 + 5X% =0 (-a_). (4.59) 



at all points with r = |(cr + + cr_) > 0. 
Since 

-l 



and 



a t du + = 8/iL-^ln |Wi| - Ci(27r|Wi|)}, (4.60) 

a L du+ = 8/iL" 1 [{Si(47r?i + ) - Si(27rD + )}cos47TM + 

-{Ci(47ra + ) - Ci(27r£> + )} sin 4™+ - 2icWi] (4-61) 

J a N du+ = 8 y uL- 1 [{Si(47rn + ) - Si(27rD + )} sin47TM + 

+{Ci(47TM + ) - Ci(27rD + )} cos4vtm + - ln(47TM + ) + lnD+], (4.62) 

it is possible to write down explicit expressions for 8X1 and 5X1 at all points on 
the perturbed trajectory. These in turn can be used to verify directly that the 
gauge constraints (|4.54|) are satisfied at all ordinary points on the trajectory. 
For example, 

\r] ab X a _{aJ) f + a b (9,a^d9 = \ f + a l da + - \ f + a N da + 

J —CT— J —(7— J — (7 — 

= -\{h cd X<LX d _)\ T=Q {a + ) + \{h cd X1X d _)\ T=0 {-a_) 

(4.63) 
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and 

r/ a6 X«((7_)^| r= o(- f 7-) = -l(h cd X c _X d )\ T=0 (-a_) (4.64) 
and so rj ab X a :((T_)5X a : (cr + , er_) = — |(/i c dX£X^)(<7 + ) as required. Similarly, 

f aV+, 0) d9 = \r[a\a + ) - sgn(<r + )a> + )] = (4.65) 

J —<7-\- 

and 

Vab X a + (a + )5X a + \ T=0 (a + ) = -l(h cd X c + X d + )\ T=0 (a + ) (4.66) 

and so r ?afe X^( ( T + )5X«( C r + ,a_) = -i(/i c ,X;X^)( ( r + ). 

To track the cumulative change in the trajectory over a single oscillation period 
At = t p = L/2, all that is necessary is to calculate the position X a + 5X a of 
the loop and the tangent vectors X± + 5X± on the spacelike section r = L/2, 
turn off the metric perturbation h a b, then evolve the trajectory forwards in flat 
space subject to the new initial conditions. The section r = L/2 corresponds to 
a± = L/2 ± a, or equivalently (since a a and 5X± are periodic functions of their 
arguments with period L) to er± = ±(er — L/2). So the perturbed null tangent 
vectors after an oscillation period are 

r L/2 

5X a + \ T=L/2 (a) = \ / a a (a-L/2,6)d6 + 5Xl\ T=0 (a-L/2) (4.67) 

J -L/2 

and 

rL/2 

8X a _\ T=L/2 (a) = \ a a (6,L/2-a)d9 + 5X a _\ T=0 (a- L/2). (4.68) 

J -L/2 

The only dependence of a a on cr_ is through the explicit trigonometric func- 
tions in (|4.48J) . and so 

r L/2 

\ a a {a - L/2,9)d9 =\L{a\a z i}\ (4.69) 

J-L/2 

where a 1 and a z are now understood to be functions of u + = (a — L/2)/L. Also, 
in view of (tOTCjl and (gS2D, 

i-L/2 /-1/2 

i / a t da + = \L I a t du + = -4/i[ 7E + ln2vr - Ci(27r)] « -9.756//, (4.70) 

J-L/2 JO 
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rL/2 

\ I a L da + = -4/x[2tt - Si(2?r)] « -19.460// (4.71) 

J-L/2 

and 

/•L/2 

| / a N da + = -4/i[ 7E + \n2n - Ci(2n)] w -9.756//, (4.72) 

J-L/2 

where 7 E « 0.5772 is Euler's constant as before. By inspection, \ f^u 2 °- z dcr + = 
0. 

Hence, if AXj_(a) = 5X±\ T=L / 2 (a) — 5X±\ T=0 (a — L/2) are the net changes in 
the perturbations 5X± over an oscillation period then 

AX» = \L[a\u + ), a z {u + ) zf (4.73) 

and 

AX!(it) = [— k, (— Acos27r-u_ + /tsin27ru_)x — (Asin27r-u_ + /tcos27ru_)y] a , 

(4.74) 

where u + = (a — L/2)/ L in the first equation and w_ = —(a — L/2)/ L in the 
second, while k = 4/x[7 + \n2ir - Ci(27r)] and A = A/j,[2n - Si(27r)]. 

The perturbed solution can now be generated by integrating the Nambu-Goto 
equation forward from the spacelike slice r = L/2 with the new initial data 
(X* + 8X+)\ T=L / 2 (o) and (Xt + 5Xt)\ T=L / 2 {a), and h ab assumed to be zero. 
Since <5X+| T=0 = 5Xt\ T=0 , the tangent vector to the initial slice is: 

T» = (x: + 5X:)\ T=L/2 = (Xl + 5X a + )\ T=L/2 -(X a _ + 5X a _)\ T=L/2 

= ±[0,sin2vra_x- cos27ra_y + sgn|w+|z] a + AX^ - AX^ (4.75) 

The vector in the span of (X£ + 5X1) \ T=L/2 and (X± + 5X1) \ T=L/2 that is orthog- 
onal to (X" + 8X%)\ T=L / 2 (to linear order in 5X±) is: 

N a (a) = (X*+8X$)\ T=L/2 + (X a _ + 5X»)\ T=L/2 

-8(X_ • 5X^X a + )\ T=L/2 - 8(X+ • 5X + Xt)\ T=L/2 
= \ [2, — sin 27ra_ x + cos 2nu^ y + sgn| u+ \ z] a 

+AX a + - 8(X + ■ AX + X°:)\ T=L/2 + AX" - 8(X_ • AX_ X a + )\ T=L/2 
+5Xl\ T=0 (Lu + ) - 8(X a _X + )\ T=L/2 ■ 5X + \ T=0 (Lu + ) 
+5Xl\ T=0 (Lu + ) - 8(X^X„)| r=L/2 • 5X_\ T=0 (Lu + ) (4.76) 
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Note here from (|4~57|) that X±\ T=L/2 ■ 5X±\ T=0 (Lu + ) = -\(h cd XlX^)\ T=0 (Lu + ) 
and so the terms involving <5X"| T= o and 5X®\ T= o in N a all cancel, leaving 

N a (a) = i[2,-sin27ra_x + cos27ra_y + sgn|u+|z] a + AX^ + AX a (4.77) 

as X + ■ AX + = X- ■ AX_ = on the slice r = L/2. 

Furthermore, the net displacement of the loop over an oscillation period can 
be measured by tracking the motion of the centre-of-mass of the loop, which for 
any constant-r slice is defined to be 

,L/2 

x (r)=L- 1 / [X(r, a) + 5X(r, a)] da. (4.78) 

J -L/2 

In view of (|4.3jh and the fact that (h ab X+X+)\ T= Q and (/i afe X"X^)| r=0 are peri- 
odic functions of a with period L/2 while X + | T= o and X_| r= o are anti-periodic 
functions with the same period (and so X r | r=0 and 5X r | T=0 are also anti-periodic), 

x (0) = §Lz and x[,(0) = 0. (4.79) 

The subsequent motion of the centre-of-mass is governed by the equation 

r L/2 

^{t) = L- 1 / (X + 5X), TT da 

J -L/2 

f L/2 

= LT 1 \ [(a L cos 27r-u_ — a N sin 27tm_)x + (a L sin27ru_ + a N cos 27r-u_)y + a z z] da 

J -L/2 

(4.80) 

where the second line follows because (X a + 5X a ), TT = a a + (X a + 5X a ) )(Ta and 
X CT + 5X CT is a periodic function of a on [—L/2, L/2}. For fixed values of r the 
acceleration components a L and a N are periodic functions of a with period L/2, 
whereas a z , cos27tm_ and sin27ru_ are anti-periodic with period L/2. So Xg(r) = 
0, and the centre-of-mass remains unperturbed by the back-reaction, a feature 
which is in any case obvious from the symmetry of the ACO loop. 

Let X a (r, a) now denote the flat-space solution generated by the initial condi- 
tions = T a (a) and X / = N a (a) on the spacelike slice r = L/2. The equation 
of motion X a , TT = X a , aa is then satisfied by taking 

X a + = 1 1 [N a (a + )+T a (a + )] (4.81) 
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and 

X a = i^ N a(_ a _^ _ T a (-a_)] (4.82) 

with a± = t ± a as usual. The position vector X itself can be found on each 
constant-r section by integrating the equation X CT = X + — X_ subject to the 
constraint x (r) = \Li. 

However, there is no simple relationship between the gauge coordinate r in 
the new solution and the Minkowski time coordinate t = X°, as 

XX=l,2 = N°(a) = 1 + \La\u + ) - k £ 1. (4.83) 

In order to facilitate comparison between successive stages in the evaporation of 
the loop, it is useful to realign the coordinates a + and cr_ by replacing them with 
two new coordinates 

ct+ = f + [N°(a) + T°(a)} da = a + + \L [ + a\0/L) d9 (4.84) 
Jo Jo 

and 



r° 

° =- lN°(a)-T°(a)}da = (l-2K)a_, (4.85) 

J — (7 — 

defined so that dX t /da + = dX t /da + = | at all points on the new trajectory. 
Note that because a± is a function of a± only, the form of the flat-space equation 
of motion X a , + _ = is preserved. 

Furthermore, since X± and AX± are periodic functions of a± with period L, 
the net change in a + or cr_ over a parametric period is: 

Aa ± = L + \ \ \ a f da+ da- = L — 2kL. (4.86) 
Jo Jo 

Thus, in a gauge coordinate system aligned with the Minkowski time coordinate 
t, the parametric period changes over a single oscillation period t p from L to 
V = L + AL, where 

AL = -2kL w -19.501/iL. (4.87) 

Since the total energy of the string loop is E = fiL, the energy radiated by the 
loop over a complete parametric period 2t p is AE = 2fiAL m — 39.002/i 2 L, in 
agreement with the original calculation of Allen, Casper and Ottewill [TH] . 
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In terms of the new gauge coordinates a + and er_, the initial data on the 
constant-time slice r = L/2 becomes 

dX«/da + = (Xl d ^)\ T=L/2 
= {§[1, sgn(a - L/2) z} a + \L[a\u + ), a z {u + ) i] a }{\ + iLa'K)}- 1 
= §[l,sgn(<7-L/2)z] a (4.88) 

as a z {u + ) = sgn(w + )a i (-u + ); and 

dX a /da^ = i[l,{-2A(l-2/t)" 1 cos27TM_ -sin27TM_}x 

+{-2A(l - 2k)" 1 sin 2ttw_ + cos 2vra_)}y] a . (4.89) 

The final step in demonstrating that the ACO loop evaporates by self-similar 
shrinkage is to show that the tangent vectors ()4.88|) and ()4.89|) are the same 
functions (to order fi) of a + and a_, respectively, as X? and X" are of o in the 
initial configuration ()4.52|) and ()4.53|) . save for a change in parametric period from 
L to L', and a rotational phase shift in (J4.89)) . 

In the case of (J4.88|) note that sgn(cr — L/2) is the restriction to the interval 
(0, L) of its L-periodic extension. As a varies from to L on the constant-time 
slice r = L/2, a + varies from L/2 to 3L/2 and a + varies from L/2 — nL to 
3L/2 — 3kL. Thus, the periodic extension of 

sgn((j - L/2) = sgn(a+ - L) (4.90) 

with period L is the same as the periodic extension of 

sgn(a+ -L + 2kL) = sgn(a+ - V) (4.91) 

on {L/2 — kL, 3L/2 — 3kL) with period V — (1 — 2k)L. But this in turn is the 
same as the //-periodic extension of sgn((X + ) on (— Z//2, 1//2), and 

«9X7&x + = ±[l,sgn(a + )zr, (4.92) 

as required. 

Referring now to (|4.89|) . the constants k and A are both of order fi, while 

n_ = = a_L-\l - 2k)- 1 = b_/L' (4.93) 

on the slice t = L/2. So 

dX a /da_ w |[1, -(2Acos2ttu_ + sin 2ttm_)x + (-2A sin 2ttu_ + cos27r?'i_)y] a 

(4.94) 
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where tt_ = a_/L'. Note here that 

r ]ab dX a /da^dX b /da- = -A 2 (4.95) 
and therefore dX a /da_ remains a null vector to first order in \i. In fact, 

- (2A cos 27ra_ + sin 2itu-) = -A sin(27™_ + ip) (4.96) 

and 

- 2A sin 2vra_ + cos 2vra_ = A cos(2vra_ + ip) (4.97) 

where 

ip = tan -1 (2A) sa 2A and A = (1 + 4A 2 ) 1/2 « 1 (4.98) 
to first order in //. It follows that 

dX a /da_ w | [1, - sin(27TM_ + 2A)x + cos(27ra_ + 2A)y] a (4.99) 

and dX a /da_ is the same function of u_ as X a ,_ is of cr/L in (|4.53|) . apart from 
the phase shift of 2A w 38.92^. 

Note finally that if ()4.88|) and ()4.89|) are used as initial data for a new oscillation 
of the loop, with f rezeroed so that the initial surface t = L/2 corresponds to 
t = 0, then 

dX a /da\ f=0 = {dX a /d& + - dX a /da^)\r = o 
= \ [1, sgn(a) z] a - i [1, - sin(-27^a/L , + 2A)x + cos(-27r<r/L / + 2A)y] a 

(4.100) 

and so 

X a (a)\ i=0 = [0,(i|a| + ||AL|)z] a +g[0,cos(-27ra/L'+2A)x+sin(-27ra/L / +2A)y] a 

(4.101) 

Hence, the kink point at a = corresponds to the spacetime point 

X a = [L/2, ||AL| z} a + £[0, cos(2A)x + sin(2A)y] a (4.102) 

which is in the same spatial position as the kink point was in the unperturbed 
trajectory at the slightly later time r = —o = XL/2n (except for a uniform 
contraction towards the centre-of-mass point). 

This observation throws some light on a question that was raised in Section 
6.11 of [SI], namely whether gravitational radiation in the weak-field approxima- 
tion would be expected to force the overall pattern of a loop (and in particular 
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salient features such as kinks and cusps) to advance or precess. In the case of the 
ACO loop it can be seen that the pattern advances. The phase shift here is due 
entirely to the form of the lateral acceleration a L . If the lateral acceleration is ev- 
erywhere inwards, and the x-y position vector is a single-mode harmonic function 
of cr_ (or er+), then it is easily seen that advance of the pattern will always occur. 
But advance of the pattern does not seem to be a general feature of all loops. 

5. Conclusion 

In this paper I have generated explicit formulae for the weak-field metric pertur- 
bations induced by an Allen-Casper-Ottewill (ACO) loop of cosmic string and 
calculated the corresponding weak-field back-reaction acceleration vector. Al- 
though the acceleration vector diverges at the two kink points on the loop, it 
turns out that the net acceleration and radiative energy loss of the loop are finite. 
Using a method first described by Quashnock and Spergel, I have determined the 
net effect of the back-reaction on the motion of the loop over a single oscillation 
period, and shown that the new, perturbed trajectory of the loop is identical to 
the original, save for a uniform contraction of scale and a small rotational phase 
shift. The ACO loop, which is rigidly rotating and has the lowest radiative ef- 
ficiency of any known loop solution, therefore evolves by self-similar evaporation 
at the weak-field level. This in turn means that, in any ensemble of loops with a 
given energy E, the ACO loops (and presumably any loops similar in shape) will 
be the most stable and longest-lived. 

The conclusions drawn here depend of course on the assumption that the 
Quashnock-Spergel method and the weak-field approximation accurately model 
the secular evolution of the ACO loop. Although there is no reason for suspecting 
that either approximation introduces spurious features into the calculation, it 
would be a useful check to be able to generate a continuously self-similar solution 
to the back-reaction equations at the weak-field level or (even better) within the 
full framework of general relativity. 

I would like to thank the referee for some helpful suggestions. 

A. Calculation of the self- acceleration vector 

As was mentioned in Section 14.11 if the field point [i, x] is a point [r, X] on the 
string the contour T in the integral (|4.15|) for divides naturally into three 
segments, which have been labelled T~ , T° and T + . When u + > the contour T 
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consists first of T~ then T° then T + , with (v + , vJ) ranging from (— |, w_ + W\ + 1) 
to (0, U-) to to (§,«- + VFi) along the segments. On the other hand, 

if u + < then the order of the segments is first T then T~ then T + , with 
ranging from (— §,«-) to to (0, m_ + W^o) to — 1). The 

divergence in h ab at an ordinary point on the string is due to the contribution near 
(v + ,V-) = {u + ,uJ), which lies on the boundary between T and Y + if u + > 0, 
and on the boundary between T and r~ if u + < 0. 

Now, if the field point [t, x] lies off the world sheet T, the line integral for h ab 
takes the form previously given in ()4.21|) . namely 

rV rVx + 1 

/i«6(t,x) = -AfxL (A; x * a6 )| s=1 ^_ -4/iL / (A+ 1 * a6 )| s= _ 1 dv_ (A.l) 

</Vi Jv 

where the limit functions Vq and V\ depend on x c = [i, x] c through equation (|4.2()|) . 
while the integrand A^ 1 ^^ depends on x c , t>_ and (in principle) v + , which is itself 
a function of x c and t>_ on V through ()4.4|) . However, the dependence of ^ ab and 
A + on v + involves a dependence on s = sgn(t> + ) only, and s is constant on the 
two sub-integrals in (jA.l|) . 

The spacetime derivatives of the metric perturbations therefore have the form 

h ab , c (i, x) = -^L{H abc {0 + , Vo) - H abc (i ^i) + H abc {-1, V + 1) - H abc (0~, V 

rV rVi+1 

+4fiL (A- 2 ^ ab A + , c )\ s=l dv. + 4fiL (A; 2 ^ a6 A+, c )| s= _ 1 ^_ (A. 

JVx Jv 

where 

A+, c = [l,-sz] c (A.3) 

and 

H abc {v + ,V) = {A-^ ab V, c )\ v _ =v . (A.4) 

Furthermore, differentiation of the light-cone condition \L{y + + V) = i — |x — 
X(Lt> + , LV)\ with v + constant gives 

^, c =2L- 1 AZ 1 [t-r,-(x-X)] c (A.5) 

with r = \L{v+ + V-) and A_ = i — r — 2(x — X) • X_ as before. 

In view of (|4.41|) . ()A.2|) and ()A.5|) . the acceleration vector a a can be written 
in the form 

a a (u + ,u.) = -8 yU L{A a (0 + ,K,)-^ a (i^i) + ^ a (-i^i + l)-^ a (0-,^o)} 

rVo rVi+1 

+8fiL / B a \ s=1 dv_ + SfiL / £°| s= _icfo_, (A.6) 

iVi Jv 
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where 

A a (v + , V) = -A-'i^X^ XZ + V a b X b _ X C + )V, C -% c X b + XI }\ v _ =v (A.7) 

and 

B a = -A~ 2 {{%X b + X c _ + %X b _ X C + )A + , C -V bc X b + X c _ V ad A + , d } (A.8) 

and again it is understood that the limit [i, x] — ► [r, X] is taken on the right-hand 
side of (CT) . 
Now, 

X\ = i[l, sz] 6 and X c _ = i[l, - sin(2vra_ )x + cos(2vrM_)y] c , (A.9) 
while the components of \]/ a & are given in (|4.17j) . Hence, 

= i(l-ss)[l,-sin(27rw_)x + cos(27ri;_)y + sz] a (A.10) 

^X b _ = i{l -cos2tt(i;_ -u_)}[l,sz] a 

+1[0, -{sin(27rz;_) - sin(27ra_)}x + {cos(2tto_) - cos(27ra_)}y] a (A.ll) 

and 

^ bc X b + XZ = 1(1 - ss){l -cos27r(u_ - it-)}. (A.12) 

Furthermore, 

A^A + , c =f and X°A+, c =±(l-ss), (A.13) 
while the general limiting values of A + and A_ at [i, x] = [r, X] are 

A+(v+,u_) = \L{u+(l - ss) + m_ - v_} (A.14) 

and 

A_(t> + , u_) = \L{u + — v + + m_ — v _ + sin27r(f_ — w_)}. (A. 15) 

For the moment, attention will be restricted to the case where u + > and so 
s — 1. It was seen in Section l4~T1 that the integral for h a b{i, x) diverges in the limit 
as [i, x] — ► [r, X], because Vo — > u_ and so A + (0 + ,Vo) — > 0, as is evident from 
(jA. 14jl . It is clear from ()A.7|) and (jA.8|) that there may be similar divergences in 
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A a (0 + , Vq) and J v ° -B a | s =i dv_. However, it turns out that both functions are in 
fact convergent, and therefore that a a is well-defined. 

To see this, note first that if u + > and v+ = + then s — s — 1 and so the 
terms (jA.lOj) and (|A.12|) vanish identically. Therefore 

A a (0 + ,V ) = -(A^ a b X b _X c + V, c )\ (0+ , Vo) (A.16) 

where, at any field point [i, x], 

X C + V )C \ {Q+ y o) = L _1 Al 1 {t — \L(v + + vJ) — s(z — |Lsf + )}|( +y ) 

= L- 1 (A: 1 A + )| (0+jVo) (A.17) 

and so 

A a (0 + ,V ) = -L-^A-JnX^ko+yo) (A.18) 

In the limit as [t, x] — > [r, X], it can be seen from (jA.15|) that (AZ )\(p+,v ) ~~ > 
2(Lu + )~ l and from (lA~TTl) that (^t X -) |(o+,v ) vanishes. Hence, A a (0+, V ) also 
vanishes in the limit. 

Similarly, because an d X+A + , c and ^/c^X^Xi are all identically zero 

when s = s, it follows that _B a | s= i is also identically zero, and therefore that 
Jy° B a \ s= i dv_ = at any field point [£, x]. Furthermore, all other limiting values 
of A + and A_ are non-zero, as V\ — > Wi(u + ) + w_ in the limit as [i, x] — > [f, X] 
and so 

A+ -> |L(2u+ - 1 - Wi) and A_ -> §L(it + - § - Wi + ^ sin 2ttWi) 
at (-§,Vi + l); 

A + — > Lu + and A_ — > \Lu + 

at (0~, Vb); and 

A+ -> and A_ -> - § - Wi + ^ sin 2^) 

at (i,V"i). It is easily verified that these limiting values are all positive if u + G 

(o,|). 

The limiting values of the three non-zero MgA b terms contributing to the 
co-rotating acceleration vector M£a b through (|A.6Jl can now be found by direct 
substitution. If u + > then 

M^A b (\, VI) -> -|L- 2 DZ 1 {(l-cos27rIVi)[l,z] a +[0,-sin27rH/ix-(l-cos27rIVi)y] a } 

(A.19) 
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M§A b (0~, V Q ) -> -iL^u-^l, y - z] a (A.20) 

and 

M b a A b (-§,^i + l) -> |L- 2 J DZ 1 (l-cos27rW 1 )[l,z] a 
-L^Dl^l, - sin2vrWZx + cos 27rW^iy - z] a 

-^L" 2 (l - cos 27rW^ 1 )D; 1 DZ 1 [0, -(1 - cos 2nW 1 )± + sin27rWiy]° 
1 

'2 

where 



+iL- 2 J D; 1 L>Z 1 W 1 [0, - sin27rW^ 1 x - (1 - cos 2nW 1 )y] a (A.21) 



D+ = 2u + - 1 - Wi and L>_ = u + - \ - W x + ± sin 2ttWi (A.22) 

The one remaining contribution to M^a b is fy^ +l MgB b \ s= _ 1 dv_. If s = — 1 
then 

%X b + X c _A + , c = i[l,-sin(27ru_)x + cos(27ru_)y- z] a (A.23) 
%X b _X c + A + , c = ±{1 -cos2tt( V _ -<)}[l,-z] a 

+\[0, -{sm(2?n;_) - sin(27ra_)}x + {cos(2tto_) - cos(27ra_)}y] a (A.24) 
% c X b + X c _ V ad A + , d = i{l - cos27r(t;_ - «_)}[!, -z] a , (A.25) 



and 
and so 



M£B b \ s= ^i = -L- 2 (v„ - 2t/L - 2z/L)~ 2 [l, -y - z] a 

+2L- 2 (v_ - 2t/L - 2z/L)- 2 [0, sin2vr(t;_ - m_)x - cos27r(t;_ - u_)y] a 

(A.26) 



Hence, 

"Vi+l 



/ 

Jv, 



M^BX^ dv. = L~ 2 {(v. - -y - z]°}|£ +1 



-2L~ 2 {(w_ -^+)" 1 [0,sin27r(t;_ -M_)x-cos27r(t;_ -w_)y] a }|^ +1 
+4ttL- 2 {[0, {Ci(27r|^_ - cos2tt(m_ - ^+) 

— Si(27r|t>_ — ip + \) sin27r(w_ — ip + )}it 
+{Ci(27r|f_ - sin27r(w_ - V+) 

- Si(27r|^_ - Wl) cos2tt( m _ - ^ + )}y] a }\v +1 (A.27) 
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where Si(x) = J* Q w _1 sin w dw and Ci(x) = — w' 1 cos w dw are again the sine 
and cosine integrals, and ip + = 2i/L + 2z/L as before. 
In the limit as [i, x] — > [r, X], therefore, 

1+1 M£B b \ s= _i dv_ - L- 2 {(2 M+ )- 1 - D- 1 }^, -y - z} a + rt^O, y] a 
+2L- 2 D- 1 [0, sin27rW^ 1 x - cos27riy iy ] a 

+4vrL~ 2 [0, (-{Ci(4vra + ) - Ci(2nD + )} cos4vra + - {Si(4vra + ) - Si(27rD + )} sin4vrM + )x 
-(Ci(4vra + ) - Ci(27i\D + )}sin47ra + + {Si(47rw+) - Si(27r£>+)} cos4vra + )y] a (A.28) 

Collecting together (lA~T9j) . (|A~20j) . (jEQ) and (jXTSSJ), plus the definition 
of W\ in terms of w+, gives finally 

"* - - z - 32ir 2 j jL- 1 F(W 1 )W 1 (l - cos 2^) (A.29) 



a 



a 



and 



a 
a 



N 



-32%fiL- 1 x 



{Si(47TM + ) - Si(27rL> + )}sin47rw + + {Ci(47ra + ) - Ci(27rL> + )} cos4vra + 
-{Si(47TM + ) - Si(27rD + )}cos47TM + + {Ci(47rw + ) - Ci(27rD + )} sin 4™ 

2nW x cos 2-irWi - sin 2irW x 
cos 2nWx - 1 + 2nW x sin 2-kW x 



+2,2ti 2 hL- 1 F{W 1 )W 1 



(A.30) 



where 



F(W) = (2n 2 W 2 + 1 - cos27riy - 2nW sin 2ttW)~\ 
Atxu + = 2tt(1 + Wi) - ^(l - cos2ttWi)/Wi 



and 



2ttD + = -1(1 -cos27rWi)/H^i. 
A similar calculation in the case where u + < (and so s 



gives: 



a 



and 



a: 
a 



-327ifiL 



-a* 



-i 



327T 2 /i^~ 1 i 71 (W / o)W / (l - cos27rW ) 



(A.31) 
(A.32) 

(A.33) 
(A.34) 



{Si(2vr(l + 2u + )) - Si(2?r.D + )} sin47ra+ + {Ci(27r(l + 2u + )) - Ci{2itD + )} cos4t™ + 
— {Si(27r(l + 2u + )) - Si(27rD + )} cos47ra + + {Ci(27r(l + 2u+)) - Ci(2nD + )} smAnu 

2tcW cos 2%Wq - sin 2kWq 
cos 2vrW / " - 1 + 2irW sin 27T1F0 



+327rVL- 1 F(W / " )W / 



(A.35) 
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where now 

Attu+ = 2irW - ^(1 - cos27tWo)/W (A.36) 

and 

2nD + = -1(1 - cos27rW )/W . (A.37) 

As can be seen from flA.29|) and (|A.34|) . a* is the same function a*(VV) of W^o 
for u + G (— 1,0) as it is of W 1 for u + G (0, |). Since = W (u + - |) for 

all it + G (0, |) the acceleration component a* is therefore a periodic function of 
u + with period |. The same is true of the lateral and normal components a L and 
a , while the vertical component a z is anti-periodic. 

B. Deriving the equation of motion from the Battye-Carter 
equation 

In this appendix I briefly describe the formalism behind the gauge-independent 
Battye-Carter equation of motion (|B.5J) . and demonstrate that the equation re- 
duces to the weak-field equations of motion (j2.13|) and (j2.14j) in the standard 
gauge. 

If (g a b,X a ) is a general solution pair to the strong-field back-reaction prob- 
lem then the projection tensor p ab corresponding to the induced 2-metric ^ AB = 
g a bX a ,A X h )B is 

p ab = j AB X a , A X b , B (B.l) 

with 7 AB the inverse of ^ab- The extrinsic curvature tensor of the world sheet T 
is therefore 

K abc = p a m p bn V n p mc (B.2) 

where V„ is the derivative operator associated with g a b. 

In terms of the derivatives of the position function X a the curvature tensor 
K abc and its trace K c = g a ^K ahc can be written as 

K abc = q^^X^A X b , B X d )CD +p^p b "T d mn ) (B.3) 

and 

K c = q c d (l CD X d , CD +p mn T d mn ) (B.4) 

with q ab = g ab — p ab the orthogonal complement of p ab , and T% c the Christoffel 
symbol associated with g a b. The equation of motion of the string, obtained by 
varying the Lagrangian density £ = — /i7 1//2 , is simply K c = 0. 
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If the equation of motion K c = is perturbed by replacing g ab with g ab + h a b 
then, to linear order in h ab , 

K c - (K abc + K a p bc )h ab + q cd V ab iy a h M - \V d h ab ) = (B.5) 

with all the geometrical quantities in this equation, except h ab , now evaluated 
using g a ° b in place of g ab [2EJ ESj • In particular, K c is no longer identically zero 
but instead of order h ab . However, the term K a p bc h ab is of second order in h ab 
and can be ignored. 

The Battye-Carter equation (jB.5|) can be specialised to the problem at hand 
by first setting gj^ = r] ab , so that = 0. Then to linear order in h ab (jB.5|) 
reads: 

q c A [n CD X\ CD -^^K.X^a X b , B X d , CD ) = -q cd l AB X a , A X b , B (h bd , a -\h ab , d ) 

(B.6) 

where q c d and 7^° are of course evaluated using rj ab rather than g ab . 

The next step in the perturbation is to replace X a with +5X a , where 5X a 
is of linear order in h ab . The only substantial change is to the term q%y X d ,cD, 
which becomes 

q c d7 CD X d 0) , CD +5q c dl CD Xf 0))CD +q c d ( 1 CD 6X d , CD -2 7 ^V D ^ 0) , A 6X b }B X d 0) , CD ) 

(B.7) 

where q c d and 7^° are now evaluated using X^ rather than X a , and 

8q c d = -5f d = -2 Vdbl AB q a b X^, A 5X a )B (B.8) 

with round brackets on spacetime indices denoting symmetrisation. 
Equation (jB.6|) therefore splits into two parts: 

q^ CD X d Q) , CD =0 (B.9) 

and 

q c dl CD 5X d , CD -2 W AB q a b X^, A 6X a , B ^ CD X d 0)}CD 
-Qdl AC 1 BD {2'U ab X^ lA 5X b , B +h ab X^ lA X b ^ B )X d ^, CD 
= —q cd 7 B X( y A X b j, B (h bd , a —^hab,d) (B.10) 

at zeroth order and first order in h ab respectively. 
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The final step is to fix the gauge. To keep the analysis as general as possible, 
suppose that the full solution pair (g a b, X a ) satisfies a gauge condition of the form 

lAB = l 1/2 KAB (B.ll) 

where kab is a constant, symmetric 2x2 matrix independent of fi, with determi- 
nant — 1. The standard gauge, for example, has kab = Vab- Replacing g a b with 
Vab + hah and X a with Xfo + 5X a then gives a zeroth-order gauge constraint 

T 1/2 lAB = KAB (B.12) 

and a first-order gauge constraint 

= Ak ab = 7~ 1/2 [2r] ab X^, ( A 5X b , B ) +h ab X( - ) , A X b (Q- ) , B 

-(o) ; c 5X b , D +lh a bX(Q^c ^(0)) 



-7 OT (^Xf 0) , c 5X fe , D +ih ab X? ),oX b D ) lAB ) (B.13) 



where, as before, all geometric quantities are evaluated using r\ ah and X? Q * in place 
of g a b and X a . 

If the zeroth-order constraint (jB.12|) is substituted into the zeroth-order equa- 
tion of motion (jB.9|) then, since q d = 8° d — p c d and 



c A COy(i -1/2/ 1/2 r7.D\ v c (T}~\A\ 

it follows that 

n „c„CD v d „,-l/2/ ■1/2„,CD v c \ CD v c (TXI^ 

as 7 1 / 2 7 C '' D = ft 073 , the inverse of kcd- The zeroth-order equation of motion is 
therefore simply 7 CD X^,cd = 0, which in the standard gauge reduces to (j2.13|) as 
required. It might seem odd that the equation q d j CD X^,cD = 0, which contains 
only two algebraically-independent components, should be equivalent to an equa- 
tion , ~f CD X^ ycD — with four independent components. However, as is evident 
from (|B.14|) . the equation r y CD X^,cD = contains information about the gauge 
choice as well as the equation of motion. 

The first-order equation of motion (jB.10|) now reduces to 

q c dl CD 8X d , CD -q c dl Ac 1 BD {2r lah Xl 0YA 5X b , B +h ttb Xfa, A X b 0) , B )Xf 0) , CD 
= —q cd/ y AB X^ lA X b ^, B (h bd , a —^h ab , d ) (B.16) 
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where 

7 AC 7 BD ( 2 ^ ab ^(o),A 5X b , B +h ab X^, A X^, B )X ( d 0) , CD 

= 7/7 7 Ak j 4 B A (0) , cd +7 (^A^jgdA ,F+2«o6^-(0))E ^-(0)^)7 X (0)'CD 

= (B.17) 
as Akab = and 7 C ' D X^,cz) = 0. So the first-order equation of motion becomes 

Qdrf CD $X d ,CD = —£[ Cd l AB X( ),A X(q),B (^6d,a — \h ab , d ). (B.18) 

Finally, after considerable algebraic rearrangement it can be seen that 

Pd7 CD 5X d , CD +p cd 1 AB X^ A X b m , B (h bd , a -\h ab , d ) 

= 7 1/ V^ 0) , b7 CD (A^c),d 

-1 AB X( 0) ,b (hat,Xfo,A +r] ab 5X a , A )^ CD X d 0))C D 

- 7 1 / 2 7 ^X ( C 0) , B j CD j EF V ab Xfo,E (X b {0) , AD Ak cf - X b 0) , FD A.kca) 

W /2 i AB x^, B l EF Xfo, E {\h ab xi, )lF + Vab sx a , F h CD (i- 1/2 i AC ), D 

= (B.19) 

(again because 7 _1//2 7ab is constant, Ak ab = and 'J CD X^, C d = 0), and so, 
given that q cd = r] cd — p cd , the first-order equation of motion (|B.18J) takes the form 

1 CD 5X C , CD = -v cd l AB X^, A X b {0) , B (h bd , a -\h ah , d ) (B.20) 

which in the standard gauge is just (|2.14|) . As with the zeroth-order equation of 
motion, the fact that (jB.20|) has four independent components does not mean that 
the system of equations is over-determined, as the tangential components (jB.19|) 
only contain information about the choice of gauge. 

Conversely, if the zeroth-order gauge condition (jB.12|) is imposed then 

k bc ^ 2 A Kab ), c = 7 1/2 ^fo)M [l CD SX c , CD +rj cd 1 AB X^ A X b Q))B (h d , a -\h ab)d )] 

W /2 (VcdSX c , A +h cd X{ 0) , A ) 7 CD Xf 0) , ra (B.21) 

and so if X? \ and 5X a satisfy the zeroth- and first-order equations of motion 
r y CD XL, C D = and (lB~20|) . the evolut ion of the first-order constraint matrix 
Ak A b is governed by the simple system of equations k bc {^ l ^ 2 Ak A b) iC = 0. 
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To analyse this system, first rotate the gauge coordinates ( A so that the sym- 
metric matrix k bc is diagonal, by setting ( A = A^'^" 4 where the orthogonal 
matrix is defined by 



A|'Ag'^ c ^^' c ' = diag(A 1 ,A 2 ). 



(B.22) 



Then the function Ua'b 1 = (Af A^ ) x ^ x ^ 2 Akab satisfies the equation k b ' c ' 
0. Furthermore, the matrix Akab is by definition (jB.13|) trace-free, and so 



/," " Ua'B'iC 



A'B' 



^k^Akab = 0. 



(B.23) 



On setting u = |(Ai£/im.' — A 2 £/ 2 / 2 /) and v = Uyy the equation k b ' c 'Ua> 
reduces to the linear system 



A'B'iC 



Au, v +Bu, 2 , = 



where 



.4 



1 
Ai 



B 



A 2 
-1 



and 



u 



(u,v) 



(B.24) 



(B.25) 



The classification and characteristics of this system of equations are deter- 
mined by the roots £ = ±a/— Ai/A 2 of the equation det(v4 — £B) = 0. Since the 
induced metric ^ AB = j~ x / 2 k ab has signature (+,—), the roots £ are real and 
distinct, and the system of equations is hyperbolic. Moreover, the characteristics 
satisfy the equation d(} / d^ 2 = £, and so have tangent directions t A ' = [£, 1]. 
Hence, 

K AB t A t B = K A >B't A 't B ' = £ 2 A^ + A^ 1 = (B.26) 

and the characteristics are null curves on the world sheet. The equation k bc {^ l l 2 Akab)-,c - 
therefore has the unique solution ^^Akab — (and satisfies the first-order 
gauge constraint (|B.13jl ) at all points on the world sheet, provided that Akab = 
on the non-null initial curve r = 0. 

It is always possible to ensure that Akab — on the initial curve by choosing 
the initial data 5X a ,A so that (jB.13|) is satisfied at all points on the curve. The 
details of this procedure as it applies to the ACO loop are explained at the begin- 
ning of Section 14.31 In the present case, the above discussion is modified slightly 
by the presence of the kinks, because j AB is undefined (that is, discontinuous) at 
the kink points, and so strictly speaking are kab and Akab (although continuous 
extensions of these functions exist trivially). Nonetheless, as is shown explicitly 
in Section I4.3( the first-order constraint Akab — is still satisfied at all ordinary 
points on the world sheet. 
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LIST OF FIGURE CAPTIONS: 

Figure 1: y-z projection of the Allen-Casper-Ottewill loop with minimum 
radiative efficiency. 

Figure 2: Backwards light cones for a field point (A) located at the kink point 
u + = U- = 0, and for a general field point (B) on the world sheet. 

Figure 3: The time component a* of the acceleration vector, in units of [i/L, 
as a function of u+. The z component a z is identical for u + > 0, but equals |a*| 
for u + < 0. 

Figure 4: The lateral component a L of the acceleration vector, in units of fi/L, 
as a function of u+. 

Figure 5: The normal component a N of the acceleration vector, in units of 
(J,/L, as a function of u + . 
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FIGURE 5 
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